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THEORETICAL PRESSURE DISTRIBUTIONS FOR SEVERAL RELATED 

NONLIFTING AIRFOILS AT HIGH SUBSONIC SPEEDS 

By John R. Spreiter, Alberta Y. Alksne, 
and B. Jeanne Hyett 


SUMMARY 


Theoretical pressure distributions on five related airfoils, 
including thin symmetrical circular-arc airfoils, in two-dimensional flows 
with high subsonic free-stream velocity are presented. The airfoils have 
various locations for the point of maximum, thickness ranging from 30 - to 
70 -percent chord and are of arbitrary, although small, thickness ratio. 

The results are obtained by approximate solution, through an iteration 
process, of a nonlinear integral equation derived from the equations of 
transonic flow theory. It is shown that the pressure distributions dis- 
play most of the principal phenomena observed in experimental studies, 
and are in good correspondence with those calculated by other methods for 
subcritical Mach numbers and for Mach numbers near 1. 


INTRODUCTION 


The equations governing transonic flows are known and well established 
by favorable comparisons with experiment. Although the difficulties aris- 
ing as a result of the nonlinearity and mixed character of the differential 
equation for the potential have hindered the advancement of the analysis, 
approximate methods are gradually emerging that permit the theoretical 
prediction of pressure distributions on a wide variety of shapes of 
aerodynamic interest. 

One of these methods is that described in reference 1 in which the 
differential equation of transonic flow theory is recast into the form 
of a nonlinear integral equation and approximate solutions are sought by 
application of an iteration procedure. The iteration procedure is of the 
successive approximation type, but differs from the related procedures 
customarily employed to determine higher approximations to the solutions 
of problems of compressible flow theory in that the quadratic nature of 
the integral equation is recognized aid retained throughout the analysis. 
This method is described in general in reference 1 and is applied to a 
nonlifting symmetrical circular-arc airfoil for which pressure distribu- 
tions are calculated for a range of Mach numbers extending from well below 
the critical. Mach number up to unity. The results indicate most of the 
principal phenomena observed in experimental studies. Certain aspects of 
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the results for free-stream Mach numbers near unity are not as satisfactory 
as for lower, hut still supercritical, Mach numbers, however, and it is 
now evident that one of the simplifying approximations needs refinement 
before results of good quality can be calculated for flows in which the 
free-stream Mach number approaches unity. 

Although it would be interesting and worthwhile to develop the 
necessary refinement to extend the Mach number range over which the method 
of reference 1 can be applied satisfactorily, the need for such a refine- 
ment has been greatly reduced by the development of a second approximate 
method described in reference 2. The latter method enables the determina- 
tion of expressions in closed analytic form for the pressure distributions 
on a wide variety of airfoils in flows with free-stream Mach numbers equal 
to or near unity. A large number of specific results are given in refer- 
ence 2, but particular attention is focused on a family of airfoils, 
including symmetrical circular-arc airfoils, that have the point of max- 
imum thickness located at 30 - , 40-, 50- , 60-,' and 70-percent chord, and 
for which experimental pressure distributions are available from refer- 
ences 3 and 4. These results indicate pronounced effects of airfoil 
shape, but provide no information on the transition from the pressure 
distributions of subsonic type that occur at all Mach numbers below the 
critical to the pressure distributions of sonic type that occur at free- 
stream Mach numbers near unity. Theoretical pressure distributions can 
be calculated for most of the intervening Mach number range by applica- 
tion of the method of reference 1, however, and it is the principal pur- 
pose of this paper to present and discuss the results of such calculations 
for the particular family of airfoils referred to above. A resume of the 
method of reference 1 precedes the presentation of the results. 


PRINCIPAL SYMBOLS 


b 
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c 
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function defined in equation (20) 
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pressure coefficient, — 
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d pressure drag 

E function defined in equation (22) 

f function defined in equations (2h) and ( 25 ) 

H - 

c 

H HtU^k-r) 173 

h half height of wind tunnel 

I f unction defined in equation (l6) 

1 ^ 2 ( 7 + 1 ) 

L 2ii L - 1 

l width of element used in approximating the chordwise velocity 

distribution 

M local Mach number 

M cr critical Mach number 

Mpo free- stream Mach number 

p static pressure 

p^ free-stream static pressure 

r 2 n/(x-I) 2 + (z-T) 2 


t maximum thickness of profile 

U TO free-stream velocity 

u,w perturbation velocity components parallel to the x,z axes, 
respectively 
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Cartesian coordinates 



NACA TN 4l48 


k 

x 

"z 

Z 

Z 


P 

r 

AOp 


6,5 


6 


T 


T 

9 

9 


a 

b 


x 

pz 

ordinates of -wing profiles 
kUoo z 
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n / i ^? 

ratio of specific heats, 1.4 for air 
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variables of integration corresponding to x,"z 
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free- stream density of air 

t 

c 

P 3 T 

perturbation velocity potential 



Subscripts 

values ahead of shock wave 
values behind shock wave 
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er conditions associated with, the critical Mach number 
L values given by linear compressible-flow theory 
W values at the wing surface 

VJT wind tunnel 

Superscript 

* values associated with the sonic point 

DESCRIPTION OF AIRFOILS 


The airfoils considered in this paper are exactly those for which 
experimental data are given by Michel, Marchaud, and Le Gallo in refer- 
ences 3 and 4. The ordinates Z of the airfoils that have the point of 
maximum thickness at 60 - and 70 -percent chord are given by 


Z 

c* 




( 1 ) 


where n has the values 3-38 and 6.05, respectively, and A is related 
to the thickness ratio t by 


n n/(n-i) 

2(n-l) 


( 2 ) 


The or di nates of the airfoils that have the point of maximum thickness 
at 30 - and 40-percent chord are given by 


Z 




( 3 ) 


where n has the values 6.05 and 3*38, respectively, and A is related 
to the thickness ratio by equation (2) . The ordinates of the airfoils 
that have the point of maximum thickness at 50 -percent chord are special 
cases of either of the families of airfoils defined by equation (l) or 
( 3 ) when n is equated to 2, and are the same as those of the airfoils 
for which theoretical results are given in reference 1. The latter 
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airfoils may be referred to as either parabolic-arc or circular-arc 
airfoils since the two classes of shapes are indistinguishable to the 
order of accuracy of the small-disturbance theory of transonic flow. 

RESUME OF THEORY 


The new results presented in this paper are determined by direct 
application of the method of reference 1 to the specific airfoils 
described in the preceding section. A resume of the parts of reference 1 
that pertain specifically to the determination of pressure distributions 
on symmetrical nonlifting airfoils is given here in the interest of com- 
pleteness, and in order to provide an opportunity for additional comments 
on certain points beyond those provided originally in reference 1. 


Fundamental Equations and Boundary Conditions 


Consider the steady flow of an inviscid compressible gas past an 
arbitrary symmetrical nonlifting airfoil of chord c and thickness 

ratio t. Introduce Cartesian coordinates x 
2 and z with the x axis parallel to the direc- 

tion of the free stream as illustrated in 
sketch (a) . Let the free-stream velocity and 
' x density be and p^, the perturbation poten-. 
tial be cp, and the perturbation velocity com- 
ponents parallel to the x and z axes be cp x , 
or u, and cp z , or w, where the subscript 
indicates differentiation. The boundary condi- 
tions require that the perturbation velocities vanish at infinity, and 
that the flow is tangential to the wing surface. The former condition 
indicates that cp is a constant at infinity. The latter condition can 
be approximated for thin wings by 



Sketch (a) 


(<P Z ) 


= u B 

z=o 00 dx 


w 


where Z represents the ordinates of the airfoil upper surface. The 
differential equation for cp in the small-disturbance theory of transonic 
flow is 


(l-^ 2 )cp xx +cp zz 


2 7+1 

= — q’x'Pxx = kc M> 

v->00 


X+XX 


( 5 ) 
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where is the Mach number of the undisturbed flow and y is the 

ratio of specific heats (lA for air). 

Equation (5) is, of course, valid only in regions where the necessary 
derivatives exist arid are continuous. Since these conditions do not hold 
where shock waves occur, and since shock waves are a prominent feature of 
most transonic flows, an additional equation is needed for the transition 
through the shock. The f undamen tal properties of a shock surface require 
that the normal component of velocity be discontinuous and the tangential 
component, and therefore cp, be continuous. The necessary relation 
follows from the classical expression for the shock polar which, in the 
small-disturbance transonic theory, is approximated by 


(l-^Cua-u^^Cwa-Wfa) 2 = M*, 2 Z±i (u a -n- b ) Z (6) 

u oo 

where the subscripts a and b refer to conditions ahead of and behind 
the shock. 

Upon solving the above boundary-value problem for the perturbation 
potential, one may determine the pressure coefficient by means of the 
approximate relation 


Cp 


P“Poo 



2_ ckp 

Uoo Sx 


( 7 ) 


An important quantity in the discussion of compressible flows is 
the critical pressure coefficient Cp associated with the local occur- 
rence of sonic velocity. The appropriate relation is found by combination 
of equation (7) and the expression obtained by equating to zero the total 
coefficient of cp^ in equation (5), and is 



2 ( 1-0 

07 + 1 ) 


(8) 


The analysis presented in reference 1 is based upon consideration 
of the preceding equations recast into a normalized form by introduction 
of reduced variables. The definitions of these quantities, together with 
additional useful relations, are given below. 


* 
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* -*2/3 00 


^cr " *»• 


where 


3 = */l-Hx, 2 


Moo 2 (r+1) 


U„ 


(9) 


In this way, the differential equat ion given by equation (5) reduces to 


q te HP M = % 


(10) 


or to 


(wi)%* ls - 0 ^ 

It is apparent from equation (ll) that u < 1 when the local velocity is 
subsonic, u = 1 when it is sonic, and u > 1 when it is supersonic. Sim- 
ilarly, the expression for the shock relation given by equation (6) reduces 
to 
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0%-~b) 2+ ( w a _!5? b) 2 (ua-^) 2 (12) 

If the shock wave is a normal wave and the flow is parallel to the 
x axis (i.e., w a = V 5 = 0 , "but u a ^ u^) , it can "be seen from, equa- 
tion ( 12 ) that the reduced perturbation velocity component u jumps from 
1 + A immediately ahead of the shock to 1 - A immediately behind the 
shock. The quantity u-(u 2 /2) , on the other hand, is equal on the two sides 
of the shock. The latter result is consistent with the fact that the 
quantity u-Ou 2 /2) corresponds, in the transonic approximation, to the 
mass flow, which is continuous through a normal shock. 


An Integral Equation for Transonic Flow 


Approximate solutions for the pressure distribution on symmetrical 
nonlifting airfoils in mixed or transonic flows are sought in the method 
of reference 1 by consideration of the following integral equation derived 
from the differential equations given in the preceding section by 
application of Green’ s theorem: 


TC(X,Tff) * TC L (x,z) + 


+00 

2 * JJ 2 
-00 



d| dT 


(13) 


where 

r 2 = n/(x-I) 2 +(z -^) 2 


The term that appears in equation ( 13 ) represents the values for u 

given by linearized compressible flow theory. Its values can be obtained 
by application of the following well-known expression: 


*--!!= r - 

* ° x Jo dg 



( 14 ) 


The derivation of equation ( 13 ) , presented in reference 1, requires 
the introduction of two slightly contradictory statements regarding the 
nature of the shock waves. They are: (a) All shock waves are assumed 

to lie in a plane perpendicular to the x axis, and (b) the shock waves 
are assumed to be norm a,! shock waves (i.e., normal to the local flow 
direction) . A new derivation of equation ( 13 ) has been given subsequently 
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in reference 5> however, in which it is shown that it is not necessary 
to introduce these two assumptions, and that this integral equation is, 
indeed, an exact relation within the framework of transonic small- 
disturbance theory. 

It is helpful to consider, before proceeding to the details of the 
approximate solution of equation (13) , a summary of the discussion given 
in reference 1 concerning some properties of _the nonlinear integral equa- 
tion and of its solution. The first and simplest property is that the 
solutions of equation (13) must approach those of linear theory when the 
free-stream Mach number is much less than the critical Mach number, 
since u « 1 and the terms involving the square of u become negligible 
with respect to those linear in u, thereby leaving only 




(15) 


Further discussion of the properties of equation (13) is facilitated 
by introduction of the abbreviation l/2 for the integral, so that 

c , n L + £ - | (16) 

where 


+00 



Although I is a function of u and is therefore unknown, it is 
informative to rewrite equation (13) by solving for u in terms of I 
and u L , thus 


u = 1 ± ^I-(2 E l -1) = 1 £ 'Jl-L 


(17) 


where 


L = 2u l -1 (18) 

Several points are to be observed at once with regard to equation (17). 

First of all, the discriminant must always be positive in order to 

obtain real values for u, thus ( 


I > L 


(19) 
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Furthermore, the choice of plus or minus sign determines whether the 
loc al velocities are subsonic or supersonic. A change in sign at the 
point where the radical is zero corresponds to a smooth transition through 
sonic velocity. A change in sign at a point where the radical is not zero 
results in a discontinuous Jimp in velocity and corresponds physically to 
a shock wave. S uch discontinuities are permissible when they proceed from 
supersonic to subsonic velocities (or from plus to minus sign in eq. (IT)) 
when progressing in the flow direction. Discontinuities in the reverse 
direction are inadmissible since they correspond to expansion shocks, a 
phenomenon which violates the second law of thermodynamics. 

In order to remove unnecessary comp li c at i ons and to facilitate the 
discussion, the following remarks are confined to the relations among I, 

L, u, and u L evaluated at the airfoil surface. In this way, each of the 
four functions reduces to a function of a single variable x and can be 
illustrated simply by curves rather than by surfaces. The subscript W 
is appended to each of the four variables to denote that the values are 
those at the wing surface. The Ly and Iy curves represent the compo- 
nents involved in the solution of equation ( 17 ), and Uy and up^ curves 
represent, respectively, the velocity distributions given by transonic 
theory and by linear theory. 

The linear theory solution up^ for subsonic flow around a nonlift- 
ing airfoil of given shape can be readily derived through application of 
equation (l4) . The corresponding values for the Ly curve follow directly 
from equation (18) . Wot very much can be stated at this point about the 
values for Iy, except that they depend on the distribution as well as 
magnitude of u and that the inequality mentioned above must be satisfied. 
The relation between the % and Ly curves is of utmost importance, 
however, and will be discussed qualitatively in the following paragraphs. 


It follows from the discussion presented in the preceding paragraphs 
that the qualitative features of the relations among the Iy, Ly, uy, 
and up^ curves are of essentially one form’ for flows in which the £ree- 

strearn Mach number is less than the critical, and of another form for flows 
in which it is greater than the critical. If the Mach number is less than 
the critical, the Iy and Iy curves never touch. Sketches of typical 
curves for this condition are shown in 
sketch (b) . The amplitudes of all 
four curves increase as the free- 
stream Mach number increases, however, 
and the Iy and Ly curves finally 
touch at one point as the free-stream 
Mach number becomes equal to the 
critical Mach number. If the free- 

stream Mach number is greater than M C o < M a0e , 



the critical, the Iy and Ly curves 
must continue to touch at at least 


Sketch (b) 


one point along the chord in order to avoid the occurrence of forbidden 
expansion shock waves. The two curves may also be tangent at a second 
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point along the chord, in "which case 
the flow may decelerate smoothly from 
supersonic to subsonic velocities. All 
present indications are that such flows 
are very exceptional if they exist at 
fl.~l.l- , and that, in most cases, the flow 
decelerates abruptly from supersonic to 
subsonic velocities through a shock wave. 
For these cases, the Iy and Ly curves 
are tangent at only one point along the 
chord and the four curves appear quali- 
tatively as illustrated in sketch (c). 
position of the sonic point and the 
point at which the sign in equation (17) changes from minus to plus. At 
some point downstream of the sonic point, the flow decelerates abruptly 
to subsonic velocities and the sign in equation (17) changes back to minus. 
If the Iy and Ly curves are continuous at this point, the Jump in uy 
is from 1 + A to 1 - A, corresponding to that of a normal shock wave, 
and (d%/dx) a = -(duy/dx)^. 

Simplification of the Integral Equation 



Sketch (c) 

The point of tangency locates the 


Since no general methods are available for the solution of integral 
equations, such as equation (13), that are both nonlinear and singular, 
recourse must be had to approximate methods. Approximate solutions of 
this equation could conceivably be worked out numerically by starting 
with a two-dimensional grid of suitably selected values for u(|,^") and 
iterating until convergence is obtained. Such calculations would pro- 
ceed by inserting the assumed values for u into the double integral 
and solving to obtain the next approximation for u(x,z) , mak in g use of 
tangency condition on the surfaces or functions represented by I and L 
as discussed in the preceding paragraphs . If the first approximation 
for u is taken to be u^, it Is shown in the appendix of reference 1 
that the final result for u, in the limit of an infinite number of iter- 
ation steps, coincides with that provided by classical iteration methods, 
provided the free-stream Mach number is less than the critical. It is 
shown, furthermore, that the power series expansions for the higher 
approximations to the solutions of the equations of transonic flow theory 
obtained by application of classical Iteration methods diverge if the 
free-stream Mach number is greater than the critical. 

The method of reference 1 does not seek approximate solutions for 
supercritical flows by taking u L to be the first approximation, but 
follows an idea of Oswatitsch expressed in references 6 and 7 that approx- 
imate solutions for mixed flows containing shock waves can be obtained if 
the starting u distribution contains discontinuities compatible with 
the shock relations. This idea is combined in reference 1 with an itera- 
tion procedure that Incorporates the considerations described above regard- 
ing the relation between the I and L functions. It is shown further in 
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reference 1 that this procedure converges rapidly and that it is not 
necessary to he highly accurate In the selection of initial values for u, 
provided care is exercised in fulfilling the tangency requirement on the 
I and L functions. 

A source of complication in the numerical solution of equation ( 13 ) 
hy an iteration process is the double integral. It is assumed in the 
method of reference 1 that approximate knowledge of the velocity distri- 
bution in space is sufficient for dete rmin ing a wor king approximation for 
the double integral. The possibility of introducing simplifying approxi- 
mations in an integral, with the effect that the error in the result may 
be rather small, appears to be the principal advantage of approaching the 
nonlinear problems of compressible flow theory, and perhaps other problems 
as well, through consideration of an integral equation rather than through 
direct consideration of a differential equation. 

The particular assumption introduced at this point in the method of 
reference 1 is that a sufficiently good approximation to the velocities 
in the flow field surrounding the airfoil can be expressed in terms of 
the local coordinate "z, the ordinates of the airfoil surface Z(Y) , 
and the desired, but unknown velocity distribution uw(x) on the airfoil 
surface. This permits one integration to be performed, thereby reducing 
the double integral of equation (13) to a single integral. The particular 
expression used in reference 1 and in the calculations of additional 
results reported in this paper is one suggested originally by Oswatitsch 
in references 6 and 7 * It possesses the properties that u starts from 
the value u^ at the airfoil surface (l = 0) with an initial rate of 
change given by the irrotationality condition and vanishes at great 
distances as l/z 3 , and is explicitly the following: 

U(X,2) = 

zfb) f 

where 

* Jb!_ 

Substitution of equation ( 20 ) in the double integral of equation (13) 
permits integration with respect to Thus, by performing this integra- 

tion and setting 1=0, the following approximate integral equation is 
obtained for u^: 



( 20 ) 
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The function E is 
E 


(^) = E(X) | jxK^-lOX^X 4 ) - 

\ b / jr(l+X 2 ) 5 L ^ 


(l-lOX^X 4 ) Zn|x| - i (l+X 2 )(25-71X 2 -X 4 , -X s ) (22) 


Although the integration interval is indicated in equation (13) to extend 
from {■ = -00 to | = +co, the contribution of the regions ahead of the lead- 
ing edge and behind the trailing edge vanishes with the introduction of . 
the velocity profile given by equation (20) since b is Infinite. The 
integration need, therefore, be carried out only over the chord c. 

Attention is called to the fact that the approximate relation for 
u(x,z) given by equation (20) is not entirely satisfactory. Evidence of 
this is provided by the facts that u is indicated to be zero in the 
region ahead of the leading edge and behind the trailing edge where b 
is infinite and that the discontinuities in u at the shock waves are 
consistent with the shock relations only at the surface of the airfoil. 

A further shortcoming of equation (20) is that it cannot be used in a 
region where b is negative since it leads to a spurious infinity when 
inserted into the integral of equation (13). For the convex airfoils 
considered in this paper, negative values for b are indicated by the 
above expression for small regions near the leading and trailing edges 
where uy is negative. Since estimates indicate that the contributions 
of these regions to the value of the double integral are very small over 
most of the chord compared with those of the regions in the middle of the 
airfoil, the contributions of regions where u is negative have been 
consistently disregarded in the present calculations. 1 

Negative values for b also result in regions of positive u^ if 
the adjacent surface of the airfoil is concave rather than convex. This 
situation does not arise in the cases considered in this paper because 
all of the airfoils have convex surfaces, but is likely to occur in other 
cases, a typical example being that of supercritical flow past an airfoil 
with a cusped trailing edge. In this particular case, the positive values 
of uvj and the chordwise extent of the region of negative b may be 
large, and the simple expedient of disregarding the influence of such 
regions could hardly be expected to succeed. - It would appear necessary, 
before pressure distributions are calculated for such cases, to reconsider 
the question of a suitable velocity profile and to introduce some other 


1 This point was improperly observed in the calculation of the results 
for the symmetrical circular-arc airfoils given originally in reference 1. 
The results for this airfoil that are presented herein have been corrected 
to conform with the above statement, even though the difference between 
the two sets of results is generally quite small . 
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expression more appropriate than that given "by equation (20) . The same 
conclusion applies also to airfoils with extensive flat regions for which 
b is infinite. 

The error 81 incurred in the value of I by the introduction of 
the relation between u and z" given by equation (20) is , of course, 
equal to the difference between the exact and approximate values for I, 
that is. 


+00 +CO +CO 

SI = Iex“-^-app lex* 7 “j/^lapp 0 " ^ ( ( lex” c lapp) tJ (23) 


where 



u a 

- s 


f 

— for b > 0 

2*[1+ (T/b) ] 4 

'Sapp = \ 


0 for b < 0 


These expressions possess two properties that make it difficult to 
estimate upper bounds for the magnitudes of 1^ and I a pp. They are that 
the region of integration extends to infinity and that a is singular, 
or infinite, at | = x, £ = "Z. In actual applications, however, the first 
property does not lead to large contributions from the distant regions 
because q^ dim in ishes sufficiently rapidly with distance from the air- 
foil, and lapp is chosen so as to assure the same behavior. The second 
property remains a factor, however, because neither q^ y nor q^pp vanishes 
at the singular point of o’. Although it is one of the principal proper- 
ties of the integral equation given by equation ( 13 ) that serious diffi- 
culties arising from this source are averted by the fact that the values 
of o change from positive to negative in a cyclic manner as the point 
| ,X moves around the singular point x,z, it is nevertheless true that 
most of the contribution to the value of I arises from the region in the 
vicinity of x,ls. It appears plausible, however, since qapp has the 
same values and first derivatives as q^ for all points along the air- 
foil chord, except those in the immediate vicinity of the leading and 
trailing edges, that the error incurred in the value of I should be small 
for some large class of smooth convex airfoils. 


It is shown in reference 1 and verified further by examination of 
the new results presented herein that considerations based on the use of 
equation (21) together with the relation between the I and L curves 
discussed in the preceding section lead to useful results for all Mach 
numbers up to that at which the shock wave is situated at the trailing 
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edge of the airfoil. At higher Mach numbers , the approximation furnished 
by equation (21) is inadequate, particularly for the region near the trail- 
ing edge, because of the presence of a pair of oblique shock waves extend- 
ing downstream from the trailing edge. Examination of the properties of 
lex i* 1 the presence of such a shock system shows that there is a log- 
arithmic singularity at the trailing edge that exactly counteracts the 
singularity in u^, and a discontinuity in its value elsewhere along the 
shock wave. This point was not understood at the time reference 1 was 
prepared and an attempt was made there to obtain results for higher Mach 
numbers by assuming that an adequate approximation for the contribution 
of the double integral to the values of u on the airfoil could be 
attained by replacing the actual oblique shock system with a normal shock 
wave situated at the trailing edge. This procedure yielded results for 
the pressure distribution on the airfoil that appeared reasonable in spite 
of the occurrence of a logarithmic singularity in u at the trailing edge. 
This singularity arose from the fact that the approximation introduced by 
use of I a pp resulted in the loss of the singularity in I, and therefore 
provided no mechanism for the proper cancellation of the singularity 
in u L . An additional difficulty that arose in reference 1 in this Mach 
number range is associated with the fact that no direct plan for iteration 
could be devised and that approximate solutions had to be sought by using 
a procedure that might be described as partly iteration and partly trial 
and error. Although it would be desirable to develop the present approach 
so that results could be obtained for the entire Mach number range, thiB 
has not yet been done because results for Mach numbers near unity can now 
be obtained by application of the hewer and simpler approximate methods 
of reference 2. The present results, together with those given in ref- 
erence 2, cover the entire Mach number range from zero to somewhat greater 
than unity with only a small gap remaining for Mach numbers equal to or 
slightly larger than that at which the shock wave reaches the trailing 
edge. 


Numerical Evaluation of Integral 


One of the principal steps in the iteration method described in 
reference 1 is the evaluation of the integral involved in equation (21) . 

Since uy and b are generally prescribed by a 
set of numerical values rather than by analytic 
functions, a numerical technique has been used 
for the integration. The process consists of 
replacing the prescribed uy distribution by a 
stepwise approximation as indicated in sketch (d) , 
introducing a mean value for b for each of the 
rectangular elements, integrating to determine the 
contribution of a single element, and summing the 
influence of all elements for which b is posi- 
tive. The contribution of a single element of 
width l and centered at x = |^, as typified by 
the shaded area of sketch (d) , is given by 
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Upon performance of the indicated operations , the following expression 
is obtained for f: 
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where 


ll+2(3C-Ti) „ Zi-2(x-|i) 

2bl ' B 2hT" 


Thus, the integral in equation (21) is appro x! mat ed as follows: 

/¥■(¥>- X ¥-[(!©.(©] 


A graph illustrating the variation of f with 2Zj_/b^ for various 
(x-f-j) /z^ is given in reference 1 and is not repeated here. 

Most of the results given in reference 1 were calculated us ing 
only 10 elements to approximate the ulyj curve. That this number is 
sufficient for the symmetrical circular-arc airfoil was checked by 
repeating the calculations in selected cases using 20 elements and noting 
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that the results were nearly indistinguishahle on the scale to which they 
are presented in reference 1. The same general procedures were followed 
at the outset of the calculations of the pressure distributions for the 
additional airfoils included in the present investigation. It was grad- 
ually concluded, however, that the extra work incurred by the use of 20 
rather than 10 elements is largely compensated by the more rapid conver- 
gence of the iteration process, and that in many cases less total work is 
required to perform the calculations with the larger number of elements. 

As a result, all of the pressure distributions presented for the airfoils 
having maximum thickness at 30- and 70-percent chord were calculated using 
20 elements. On the other hand, the pressure distributions for the air- 
foils having maximum thickness at 40- and 60-percent chord were calculated 
using 10 elements for all Mach numbers except the critical, for which 
20 elements were used. 


Determination of uj^ 


The values for ujv. can be determined from the general solution 
for u L given in equation (l4) by performing the indicated operations 
and setting z = 0, or more simply from the following expression derived 
from equation (l4) : 


= 



dZ 

dt x-T 


(27) 


The Cauchy principal value Is used in the integration. In this way, the 
following expression is obtained for uj^ for the airfoils described by 
equations (l) and (2) : 


t n n /(n-i) 
UL W"“^T“ 




1 

n-nH-1 





+ 


( 28 ) 


where m is such an integer that 0 < (n-m) <1. If n is any integer, 
n-m is zero and the integral, disappears. Thus equation (28) reduces to 
the following simple form for the symmetrical .circular-arc airfoil for 
which n = 2: 



NACA TN 41A8 


19 


u % 




c-x 


2 


(29) 


If n is not an integer, the 
integral remains. It can he eval- 
uated analytically if n is any 
integer plus certain simple frac- 
tions. In the present investiga- 
tion, however. It is necessary to 
consider the values 3-38 and 6.05 
for n and the integrals are 
evaluated numerically. No diffi- 
culty occurs at the point where 
"f = x since the integrand is not 
singular hut is equal to 

(n-m) /(x/c) . The general 

expression for uj^ for the air- 
foils described hy equations (2) 
and (3) can also he written in a 
similar manner, hut It is suffi- 
cient to note that the results are the same, except for replacement of 
x/c hy l-(x/c), as those found hy application of equation (28). The 
chordwise variation of ur for each of the five airfoils of the present 
investigation Is shown graphically in sketch (e) . 



Iteration Solution of Integral Equation 


Approximate results for pressure distributions on thin airfoils are 
determined, in the method of reference 1, hy solving the simplified inte- 
gral equation hy use of an iteration procedure in which certain operations 
at each step are performed numerically and other operations are performed 
graphically. In particular, the value of the integral, or the function 
designated in the above discussion hy I, is determined numerically; the 
tangency condition between the I and L functions is satisfied, when 
necessary, hy use of graphical technique; and the resultant values for uy 
are determined at each step hy use of the quadratic formula as indicated 
by equation ( 17 ) . These techniques have been applied without modification 
in the calculation of the new results presented in this paper. It should 
he remarked, however, that there are two categories of problems discussed 
in reference 1 for which no new results are given in this paper. One of 
these is that of shock-free supercritical flows for which no further 
attempts were made to attain solutions because It seemed evident that the 
same negative result reported in reference 1 would again he found. The 
other category Is the one mentioned previously that includes the Mach 
number range near unity at which the shock wave is situated at the trail- 
ing edge and the approximation of the variation of u with z’ furnished 
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by equation (20) is clearly inadequate. In spite of this known deficiency 
and the lack of a direct plan of iteration, numerous calculations were 
actually made for this category. Although the quality of the results was 
generally comparable with that of the corresponding results reported in 
reference 1 for the symmetrical circular-arc airfoil, it was not equal to 
the quality of the results for lower Mach numbers. This situation, 
together with the decreased necessity for the results occasioned by the 
development of the theory of reference 2, has caused these results to be 
omitted from the present report. 

The two categories for which new results are given in this paper are 
(a) subcritical flows and (b) supercritical flows for which the shock wave 
is situated forward of the trailing edge. The details of the iteration 
method are such that the Mach number or value of g^ for a given calcu- 
lation can be selected in advance for subcritical flows, but not for 
supercritical flows. In the latter case, the position of the shock wave 
is fixed instead, and the associated Mach number is found as part of the 
solution. For this reason, results given for supercritical flows are 
selected so that the shock wave is situated at even intervals of the 
chord, rather than so that there are even intervals of Mach number or g . 


BESUITS AND DISCUSSION 
Presentation of Results 


The theory described in the preceding pages has been applied to 
calculate pressure distributions for both subcritical and supercritical 
Mach numbers on each of the five families of airfoils described at the 
beginning of this paper. As noted previously herein, the Mach number 
range investigated extends in each case from well below the critical up 
to the lowest Mach number at which the shock wave is at the trailing edge. 
The results are given in numerical form in table I and are illustrated 
graphically in figures 1 through 5* Each figure consists of two parts. 
Part (a) presents the results for subcritical flows and part (b) presents 
the results for supercritical flows. 

As indicated previously herein, the pressure distribution associated 
with the value of g for which the shock wave is situated at the trail- 
ing edge is inaccurate in the vicinity of the trailing edge. The results 
presented for this case are consequently indicated by a dashed line over 
the last few percent chord. 

In addition to the above results, each part of each figure contains 
an extra set of results indicated by a dotted line that is useful for 
reference and comparison. The dotted line in part (a) of each figure 
indicates the pressure distribution provided by linear theory for the 
lowest Mach number (or most negative g ) considered for each airfoil. 

That in part (b) of each figure indicates the theoretical pressure 
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distribution given in reference 2 for free-stream Mach number 1 (g = 0) . 
The latter results are shown in reference 2 to be in substantial agreement 
with the experimental results of references 3 4. 


Ex ami nation of the results for the various airfoils reveals that the 
pressure distributions for supercritical flows bear an interesting rela- 
tion to the pressure distributions for the critical Mach number and for 
free-stream Mach number 1. As the Mach number is increased beyond the 
critical, the shock wave increases 
in strength and moves rearward 
across the chord, but the values 
of Cp over the forward portions 
of the airfoil remain nearly the 
same as at the critical Mach numbei 
When the Mach number is increased 
sufficiently that the shock wave 
reaches the trailing edge, the 
entire pressure distribution resem- 
bles in form, but not in level, the 
pressure distribution at Mach num- 
ber 1. It can be seen, however, 
that the chordwise exbent of the 
supersonic region and the rate at 
which it grows as the Mach number 
is increased beyond the critical 
are strongly dependent on the shape 
of the airfoil. These statements 
are illustrated by the curves pro- 
vided in sketch (f) that show the 



variation with g ra of the locations of the sonic point and the shock wave 
for each of the airfoils. The dotted lines for values of g^ near zero 
refer again to the results given in reference 2 for the location of the 
sonic point at Mach number 1 and vicinity. 


The results for the airfoil having maximum thickness at 30-percent 
chord are particularly interesting because they indicate that the position 
of the shock wave changes very rapidly with small, changes in g^ or Mach 
number. It is possible that this result may have a bearing on some of the 
unsteady conditions that are observed in experimental investigations (see, 
for instance, ref. 8) . These conditions are generally associated with 
extensive shock- wave boundary-layer interaction. The circumstances 
encountered on the airfoil with the point of maximum thickness at 
30 -percent chord suggest that unsteady phenomena of major importance may 
also occur even in inviscid flows if there is a little unsteadiness in 
the main stream. The circumstances associated with this airfoil are in 
marked contrast to those for the other airfoils considered herein for 
which the position of the shock wave is more strongly dependent on g^. 

Examination of the simplified integral equation given in equation (21) 
discloses that the great sensitivity of the shock position to small changes 
in Mach number exhibited by the results for the airfoil with maximum 
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tMckness at 30-percent chord is not to "be associated, in general, with 
a forward location of the point of maximum thickness, hut with a vanish- 
ing of the curvature of the rear portion of the airfoil. This conclusion 
follows directly from the fact that the contribution to the value of the 
integral of equation (21) arising from such regions is very small because 
b is large according to equation (20), E[(g-x)/b] is large logarithmically 
according to equation (22), and the quotient E/b in the integrand van- 
ishes in the limit as the curvature approaches zero. It thus follows that 
the I curve, hence the ujuj curve, and finally the free- stream Mach 
number are very insensitive to variations in the location of the shock 
wave if the shock wave is situated adjacent to a part of the airfoil where 
the curvature is very small. Although this conclusion appears to be a 
general property of the simplified integral equation given in equation (21) 
some caution should be exercised in its immediate acceptance and widespread 
application to physical problems because the approximate relation, equa- 
tion (20), introduced in the development of equation (21) from the exact 
integral equation given in equation (13) obviously deteriorates in quality 
as the curvature vanishes. The agreement between the present results and 
those for Mach numbers near 1 given in reference 2 suggests, however, that 
the quality of the specific results given in this paper is not unduly 
impaired by this factor. 

The results given in figures 1 through 5 are presented again in __ 
figures 6 through 10 in an alternative form in which the variation of Cp 
with | OT is plotted for various stations along the chord. This form of 
presentation is the counterpart, in terms of reduced quantities, of the 
plots commonly given in reports of experimental investigations showing 
the variation of Cp with at prescribed points on the surface of a 
wing or body. An additional set of lines is included on each of these 
figures showing the theoretical results given in reference 2 for the vari- 
ation of Cp with g OT for free-stream Mach numbers near 1. It can be seen 
that these results are in essential agreement with those calculated by the 
method of reference 1 for lower Mach numbers. This agreement is of par- 
ticular interest because the two sets of results have been calculated by 
methods that are new and approximate, but that differ widely in the 
details of the approximation. 

It Is evident from examination of figures 1 through 10 that the 
results calculated by the method of reference 1 display the anticipated 
qualitative aspects of transonic flows associated with the part of the 
Mach number range to which they apply ( 1 . e . , appearance , rearward move- 
ment, and increase in strength of the shock wave as the Mach number is 
increased beyond the critical, etc.) and are in good quantitative agree- 
ment with other theoretical results for low Mach numbers and for Mach 
numbers near unity. The remaining sections of this paper are concerned 
with comparisons of. the present results with those of higher order theory 
for subcritical flows and with those of experiment for supercritical flows. 

_The results given in figures 1 through 10 are all expressed in terms 
of C p and. g^. Since the extraction of the value of that is asso- 

ciated with a given pair of values for g^ and the thickness ratio t 
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requires the solution of a cubic algebraic equation and is somewhat 
inconvenient , figure 11 has been included to give the variation of 
with for various r. It is assumed in the preparation of this graph 

that 7 is 7/5. 


Comparison With Existing Higher Approximations for 
Sub critical Flows 


Although the results for supercritical flows represent the primary- 
contribution of the present investigation, it is not without interest to 
examine the results for subcritical flows in somewhat greater detail and 
to compare them with existing results given by linear theory and by higher 
approximations. These comparisons are of particular importance in the 
present discussion because it is only for subcritical flows that a good 
evaluation of the degree of accuracy of the results can be achieved at 
the present time. Even this discussion is handicapped somewhat by the 
fact that all other theories for subsonic flow around airfoils are also 
approximate and that no exact solutions are known. 


From the present point of view, the most significant method that has 
been used to obtain higher approximations to the solution for subcritical 
flow around thin airfoils is the method of successive approximation in 
which the solution is expressed in a power series in thickness ratio. 

In this method, the first term is the result given by linear theory, and 
the coefficients of successive terms are determined by iteration. Although 
a general procedure for the evaluation of the second approximation has 
recently been given by Van Dyke 
(ref. 9 ) } the determination of 
the third approximation has been 
accomplished for only a few spe- 
cial shapes. One of these is 
the nonlifting symmetrical 
circular-arc airfoil for which 
the second approximation has 
been given by Hants sche and 
Wendt (ref. 10) and the third 
approximation by Asaka (refs. 11, 

12, and 13) . 2 Sketch (g) shows 
a comparison of the variations 
of Cp with at the midpoint 
of such an airfoil having a 
thickness ratio t of 0.10, as . 

indicated by the present calcu- Sketch (g) 

lations and by the first, second, and third approximations. It can be 
seen that the present results approach those given by the first 

2 The results for the third approximation given in the present report 
differ from those obtainable directly from the expressions given in either 
references 11 or 12 and 13 because of the correction of some misprints. 
These corrections have been verified by correspondence with the author. 
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approximation (or linear theory) for small Mach numbers, hut depart 
therefrom with increasing Mach number and are much closer to the higher 
approximations for Mach numbers near the critical. At Mach numbers 
greater than the critical, the results of the present calculations differ 
markedly from those indicated by the method of successive approximation. 
This is proper, however, since the range of applicability of the latter 
method is • confined to subcritical Mach numbers . 

It should be noted that the curves labeled first, second, and third 
order present the results indicated by successive approximations to the 
solution of the exact equations for inviscid compressible flow. It is of 
greater significance, from the point of view of evaluating the accuracy 
of the present calculations , to compare the present results with those 
indicated by successive approximations to the solution of the simplified 
equations of transonic small-disturbance theory. The latter results can 
be calculated by use of the following expression, which is readily derived 
from Asaka* s result by taking the limiting form consistent with the 
approximations of transonic flow theory: 


(Uocft) 
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where for the midpoint of a circular-arc airfoil 
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successive approximations . 


Sketch (h) 


In order to provide further insight into this situation, three 
additional curves are included in sketch (h) to show the variations 
indicated "by the following expression: 


auo 1/3 c 

t 2/3 % 


= % - 2! » 




where for the midpoint of a circular-arc airfoil the coefficient, a~, is 
given in equation (31) , Notice that 



The foregoing expressions represent the results of the Nth approximation 
obtained by application of an alternative method of successive approxima- 
tions that involves the solution of quadratic rather than linear equa- 
tions at each step of the iteration process. This method, described 
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"briefly in reference 1 and discussed at greater length in reference 2, 
yields results at each step that overestimate, rather than underestimate, 
the variation with of the peak negative values for Cp. Although 

the difficulties of integration are as great as in the classical method 
of successive approximations and only the first few steps can be evaluated 
in any specific application, it is shown in "the appendix of the report 
version of reference 1 that the results obtained for flows that are sub- 
sonic everywhere converge, in the limit of an infinite number of iteration 
steps, to the same result as ultimately obtained by application of the 
classical method of successive approximations. The result obtained by 
application of the quadratic method clearly terminates with the occurrence 
of sonic velocity somewhere in the flow field, however, and no results are 
provided for mixed or transonic flows. This termination of the result is 
in distinct contrast to the apparent behavior of the results of the classi- 
cal method of successive approximation, but closer examination of the rela- 
tion between the two sets of results reveals that the infinite series with 
which the latter results are expressed does not converge at Mach numbers 
greater than the critical and that the values indicated for transonic flows 
are false. Since the exact solution of the equations of transonic flow 
theory is presumed to indicate a variation of Cp with that is somewhat 
greater than that indicated by the third approximation of the classical 
method, but somewhat less than that indicated by the third approximation 
of the quadratic method, it is evident upon examination of sketch (h) that 
the proper trend is defined within rather narrow limits, and that the 
results of the present numerical calculations fall within these limits. 

-4.0 Although it is immediately 

evident from the foregoing dis- 
cussion that results correspond- 
ing to those illustrated in 
sketch (h) for the 10-pereent- 
“ 3 -° thick airfoil could be calculated 
0 P readily for airfoils of other 
thickness ratios, graphical pre- 
sentation of the results in 
_20 terms of Cp and would 

require a complete set of curves 
for each thickness ratio consid- 
ered. The corresponding values 
for thin airfoils of arbitrary 
-1.0 thickness ratio can be summa- 
rized concisely by a single set 
of lines, however, by plotting 
the results in terms of the 
reduced variables Cp and 1^. 

0 Such a plot is given in 

sketch (i). It can be seen that 
the relationship between the 
various approximations remains 
substantially the same as illustrated in sketch (h) for the special case 
of a 10-percent-thick airfoil. 



4.0 - 3.0 -20 - 1.0 

^a> 


Sketch (i) 
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The corresponding_re suits 
for the variation of Cp 
with | at. the station, x/c, 
along the chord at which the 
present results indicate sonic 
velocity is first attained as 
the Mach number is increased 
are shown for the other four 
airfoils of the present inves- 
tigation in sketches (j) 
and (k) . It is necessary to 
present only two plots because 
the results for suberitieal 
Mach numbers are the same for 
the airfoils having the point 
of maxim-urn thickness located 
at 30 - and 70 -percent chord, 
and similarly for those having 
maximum thickness at 40- and 
60-percent chord. Since 
approximations higher than the 
first have not been determined 
for comparison with the pres- 
ent numerical results for the 
latter four airfoils, an addi- 
tional curve is included in 
sketches (i), (j), and (k) 
illustrating the variation 
of Cp with indicated by 
the following pressure- 
correction formula derived in 
reference 2: 


Cp - 2I e 


-HM* 


(34) 



C, 

Sketch (j) 


is the value 



where Cp^ 

for Cp indicated by linear- 
ized compressible flow theory. 

The relation given by equa- 
tion (34) is approximate, but 
it can be seen from sketch (i) Sketch (k) 

that its use leads to results that are more accurate than those indicated 
by linearized theory. Although it has not been proved in general, it 
appears from examination of the results for the symmetrical circular-arc 
airfoil and other cases for which higher order approximations are avail- 
able that the variation of Cp with indicated by equation (34) is 
somewhat greater than the proper trend. 
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It can be seen from the foregoing comparisons that the present 
numerical results appear in each case to be in good accordance with other 
theoretical results for subcritical Mach numbers. There are no correspond- 
ing theoretical results available with which the present results -for super- 
critical Mach numbers can be compared, however, and evaluation of the 
accuracy of the results in this range can only be attempted at the present 
time by comparison with experimental results. 


Comparison With Existing Experimental Results 
for Supercritical Flows 


The appeal to experiment to judge the quality of an approximate 
solution of the equations of compressible flow is nowhere a more uncer- 
tain procedure than in the range of Mach numbers somewhat greater than the 
critical in which the shock wave is situated upstream of the trailing edge. 
The principal reason is that the theoretical results are based on the 
assumption of an inviscid gas, whereas numerous experimental investigations 
(see, e.g., refs. 14 through 19) indicate that viscous effects associated 
with interaction of the shock wave and the boundary layer are frequently, 
if not invariably, of substantial magnitude in the general vicinity of 
the point where the shock wave meets the boundary layer. A resume of some 
of the more prominent aspects of these phenomena is given in reference 1. 


As a start toward a more specific 
discussion of the present results, con- 
sider the comparison shown in sketch ( z) 
of the theoretical and experimental 
variations of Op with g ro at the mid- 
point of symmetrical circular-arc air- 
foils. The circles joined by the solid 
line represent the calculated results 
for a thin airfoil of arbitrary thick- 
ness ratio. The other data points 
represent the experimental results from 
reference 3 for four airfoils having 
thickness ratios of 0.06, 0.08, 0.10, 
and 0.12 and should, according to the 
transonic similarity rule, define a 
single curve. Although the theoretical 
and experimental results are in reason- 
Sketch { l) able agreement at Mach numbers near 1 

and are qualitatively similar at all Mach numbers, it can be seen that 
substantial quantitative discrepancies exist for Mach numbers in the gen- 
eral vicinity of the critical. It is unlikely that these discrepancies 
at the midchord station can be attributed to shock-wave boundary- layer 
interaction. Examination of the results for the other airfoils of the 
present investigation reveals similar discrepancies when compared with the 
experimental data of Michel, Marchaud, and Le Gallo given in reference 4. 
Since the results for the circular-arc airfoil are typical of those for 
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al 1 of the airfoils considered herein, most of the following discussion 
is concerned with closer examination of the results for this particular 
airfoil. 

Experimental pressure distributions for supercritical flows past 
symmetrical circular-arc airfoils have been published not only in refer- 
ences 3 and 4, but also in references 14 and 15 by Liepmann, Ashkenas, 
and Cole, and in reference 16 by Wood 
and Gooderum. It should be noted, 
however, that the investigations 
reported in references 14 and 15 are 
concerned pr imar ily with boundary- 
layer shock-wave interaction and the 
authors include statements in refer- 
ence 14 casting doubt on the accuracy 
of the values indicated for the free- 
stream Mach number, a quantity of 
only secondary importance in their _ 
investigations. The variation of Cp 
with at the midpoint of these 

airfoils is presented in sketch (m) 
in exactly the same manner as in 
sketch ( 2) so that the two sets of " 40 _3 ° ' 2 -° " l£) 0 

results should coincide. It can be 

seen that such coincidence of the Sketch (m) 

results Is again observed at Mach numbers near 1, but that there is wide 

scatter in the results for Mach numbers in the vicinity of the critical. 

It is interesting to observe in all these cases, as well as in others for 
Mach number 1 given in reference 2, that the experimental results for the 
thicker airfoils are generally in better agreement with the calculated 
results than those for the thinner airfoils. This trend is contrary to 
what might be expected if a major source of error were associated with 
the assumption of small disturbances in the establishment of the equations 
of transonic flow theory. 

Inasmuch as the differences among the various sets of results exist 
at Mach numbers less than as well as greater than the critical. It is 
possible to explore the accuracy of the various results by comparison with 
the third-order approximation of Asaka to the solution of the complete 
equations of subsonic compressible flow theory. Since the latter theory 
includes terms disregarded in transonic flow theory and results for air- 
foils of different thickness ratios do not col lapse perfectly onto a single 
line when plotted in terms of Cp and the results are given in 
sketch (n) in terms of Cp and M^. In order to avoid confusion due to 
overlap, the intermediate thickness ratios, t = 0.08 and t = 0.10, have 
been omitted. It can be seen that none of the experimental results are 
in good agreement with the values indicated by higher order theory. Inas- 
much as there is considerable evidence from various sources that higher 
order theory generally provides results that agree well with those found 
experimentally for Mach numbers up to the critical, it would appear that 
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-•8 | ~ 1 1 n a \ 1 1 the source of the 

discrepancy should be 
sought in some facet of 
the experimental 
technique . 

Some disagreement 
between theoretical and 
experimental pressure 
distributions, although 
gener ally of much smaller 
magnitude than that 
observed in sketches ( l) 
through (n) , must always 
be anticipated to result 
from the presence of the 
boundary layer. The 
influence of this factor 
may be somewhat greater 
than usual in the expert - 

Sketch (n) mental investigations of 

references 3 and 4 because the airfoil is not mounted in the middle of the 
tunnel as is customary, but is simulated by a bump on the tunnel wall, and 
is hence imbedded in the wall boundary layer. The authors of references 3 
nnri 4 have considered this point, however, and Included some evidence that 
indicates the experimental results on the 12-percent-thick circular-arc 
airfoil at a Mach number of 0.99 are substantially the same as those on a 
similar airfoil mounted in the middle of the tunnel and having a boundary 
layer artificially made turbulent . No corresponding information is pre- 
sented for lower Mach numbers or for the thinner airfoilB for which the 
relative effects of the boundary layer would, be greater. Although it is 
recognized that the influence of the boundary layer on the pressure dis- 
tribution is not the same for all Mach numbers, it is difficult to account 
for the rtorI 1 discrepancies between theory and experiment near Mach num- 
ber 1 and the large discrepancies near the critical Mach number by 
boundary- layer effects alone. It appears, therefore, that some other , 
factor might be supplying a substantial contribution to the observed 
discrepancies . 



A source of uncertainty that is present in all wind-tunnel testing 
is associated with the finite dimensions of the test section. The study 
of wall interference and the determination of formulas for the calculation 
of corrections have consequently been the subject of numerous investiga- 
tions (see, e.g., refs. 20 through 25). It is found, for airfoils and 
test sections of typical proportions, that the influence of the walls is 
generally smal 1 at low Mach numbers, but increases substantially as the 
Mach number increases toward and into the transonic range. Preliminary 
investigation reveals, moreover, that it appears plausible to attribute a 
substantial amount of the present discrepancies to this source. 
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In order to be more specific, consider the problem of transonic flow 
past a thin nonlifting airfoil of chord c and thickness ratio r mounted 
in the center of a wind tunnel of height 2h. The aerodynamic properties 
of this airfoil are to be compared with those of the same airfoil in an 
unbounded flow with free-stream velocity U^. It is considered that the 
apparent free-stream velocity in the wind tunnel, as determined from 
measurements of pressure at a point far upstream in the test section or 
in the plenum chamber of a ventilated wind tunnel, is not necessarily 
equal to U^, but is given by + u^ where u^ is much less than U^. 

The appropriate equations for the study of such flows are thus equa- 
tions (4) through (6) together with the following relation for the 
pressure coefficient C-^ on the airfoil in the wind tunnel 



g(u-Q 


(35) 


and an additional boundary condition at z = ±h dete rmine d by the nature 
of the wa ll . This condition is w = 0 if the test section has solid 
walls, u = u^ if the test section is an open jet, and a more complicated 
relation if the test section has porous or slotted walls. It is not nec- 
essary to be more specific about the boundary condition for the latter 
case at this point, since it is shown in references 21 and 22 that the 
influence of a ventilated wall is intermediate between that of a solid 
wall and an open jet. 

The linearized theory of wind-tunnel wall interference follows 
directly from the equations enumerated above upon replac ing the right- 
hand sides of equations (5) and (6) by zero. It is found, provided the 
dimensions of the section are large compared with those of the airfoil, 
that the net effect of the walls at subsonic Mach numbers is that the 
pressure p (not Cp) and local Mach number at the surface of the air- 
foil in the wind tunnel in a flow with apparent free-stream velocity 
Uoo + Uoo is the same as at the surface of an identical airfoil in an 
unbounded flow with free-stream velocity U^. The quantity u^ is 
given by 


^oo _ £ t 

Uqo (1-L^o 2 ) s/2 H 2 


( 36 ) 


where H represents h/c and C is a constant, the value of which 
depends on the nature of the walls of the test section; for example. 
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A ' 

C = - — — solid walls 

6 tc 

> (37) 

- 

G »» + ~ — open jet 

12 tc 

— 


where A<. is the area of the airfoil section. The influence of the 
tunnel walls on Cp follows directly upon combination of equations ( 35 ) 
and ( 36 ) . It is 


- v ■ - 2 £ 

(38) 


Comparison with other results presented in this report is facilitated 
by rewriting equations ( 36 ) and ( 38 ) in terms of C p and 1^, thus 


c 

A| “ tco loo WT (t NS/Sh 2 

' 5 co' 

(39) 

« 

AC p = Gp - C^ = 2 A| ro 

(40) 


where 



8_ - - , H . (U„k T ) 1/3 H 

(u^t) 273 




If the influence of the walls is examined for a family of airfoils of 
different thickness ratio and with the apparent Mach numbers in the wind 
tunnel selected so as to maintain a given and thereby to secure 

similarity of transonic phenomena such as the critical Mach number, etc., 
it is evident from equations (39) and (40_)_ that IaI^I and |ACp| are 
inversely proportional to the square of H. It follows immediately that 
[A^ | and |ACp j diminish If either the height-chord ratio h/c of the 
tunnel or -the thickness ratio r of the airfoil is increased. 
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Added insight into the nature of 
the influence of the tunnel walls can 
he had by examining the curves given 
in sketch (o) showing the variation 
of Cp with at the midpoint of a 
symmetrical circular-arc airfoil in 
an unbounded flow (or the variation 
of wi-fc h for infinite H) 

together with the corresponding 

values for Cn with 6 computed 

°°WT 

hy use of equations (39) and (40) for 
the same airfoil mounted in either a 
closed wind tunnel or an open jet of 
such dimensions that H equals 
either 1 or -J2. The value of unity 
for H corresponds, for instance, to 
that associated with a 6 -percent- 
thick airfoil in a wind tunnel having 
a ratio of semiheight to chord h/c 
of 7/3 (the same as that for the 
experimental data shown in sketch ( l) ) 
at an effective free-stream Mach num- 
ber °f approximately O. 75 . Although the curves are continued to 

supercritical Mach numbers for which the use of linear theory to calculate 
the influence of the tunnel walls may be open to question, it is inter- 
esting to observe that the differences between the results for an 
unbounded flow and for an airfoil in an open jet are qualitatively simi- 
lar to those that can be observed in sketch ( l) between the theoretical 
results for an unbounded flow and the experimental results measured in a 
wind tunnel with slotted test section. The differences between the theo- 
retical and experimental results are larger than indicated by the linear 
theory of wall interference, but quantitative agreement should perhaps 
not be expected since it is well known that substantial nonlinear effects 
occur at Mach numbers of the order of the critical and greater. 



Sketch (o) 


An interesting property of the linear theory for the influence of 

the tunnel walls is that there is a maximum value for f , and hence 

°°WT 

the apparent free-stream Mach number, in a closed tunnel. This result 

is at least qualitatively consistent with the familiar property of choking 

of wind tunnels with solid walls. If the dimensions of the test section 

are sufficiently large compared with those of the airfoil, the values 

for Ct-, associated with the maximum value for E are on the exten- 
'■'0.™ ^WT 

- 26 ra is equal to -26, 

which is the reduced quantity associated with the local Mach number, this 
result indicates that the local Mach number distribution measured on an 
airfoil in a solid-wall tunnel of sufficiently great dimensions under 
choking conditions is the same as in an unbounded flow with free-stream 
Mach number 1. This conclusion is in complete accordance with the more 


sion of a line 




2 | = 
3 00 


"Pm _■> 


Since 


9p 
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reliable theoretical results, based on nonlinear transonic flow theory, 
for flow past a wedge airfoil and a lifting flat plate (refs. 24 and 25), 
and with experimental results for an NACA 64AD08 airfoil tested in several 
test sections of different sizes (ref. 26). 

Linear theory indicates infinite corrections when M*, = 1 and must 
be replaced by a corresponding theory based on the nonlinear equations 
of transonic flow theory. This theory is not developed as far as the 
linear theory, but it is clear from similarity considerations that the 
functional form for the expression for CL. for a family of affinely 
related airfoils in a wind tunnel is 


V - I- H) <M) 

plus additional parameters indicative of the porosity of the walls or of 
the geometry of the slots for wind tunnels with ventilated test sections. 

If it is assumed that the local Mach number distribution measured 
on an airfoil in a wind tunnel is the same as that measured in an unbounded 
flow, in accordance with the indications of linear theory, equation (40) 
again applies together with an expression for that has the following 

form: 

^co = *(!«> 5) . (te) 


Equation (39) , obtained from linear theory, is consistent with this 
functional expression and is, moreover, the precise relation that the 
solution of the nonlinear equations of transonic flow theory must approach 
as the free- stream Mach number tends toward zero. Although analysis of 
the influence of the walls based on consideration of ACp and A| ra as 
described above is similar to that successfully employed for low-speed 
flows, it may not be the best procedure to follow for Mach numbers near 1. 
The reason is that the corrections must inevitably become large at Mach 
numbers near 1 as a result of the fact that Cp - 2^ is independent 
of in this range. 

Alternative procedures that lead, to smaller corrections at Mach 
numbers near 1 are to consider ACL, at £ .= g . or to find A£ bo 

J? °°WT 

that Cp^ = Cp at some representative point on the airfoil. If the 

former point of view is adopted it follows from equation (4l) that ACp 
is a function of and H. At Mach number 1, ^ is zero and A3p 
depends only on H, thus 



=o 

00 


= f(H) 


(43) 
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If it is assumed that |ACp | decreases with increasing H, as appears 
reasonable since |ACp| would be expected to be smaller if a given airfoil 
were tested in a large wind tunnel than in a small wind tunnel, it follows 
again that {ACp| is smaller for a thick airfoil than for a thin airfoil 
of affinely related geometry tested in a given wind tunnel. 

Although the above discussion applies strictly only to families of 
affinely related airfoils, certain properties of the asymptotic solution 
for the flow at great distances from an airfoil at Mach number 1 suggest 
that the magnitude of |ACp| is greater for airfoils having maximum thick- 
ness far aft than for those having maximum thickness far forward. The 
reason Is that it is not so much the thickness and chord of the complete 
airfoil that matters in the expressions for the asymptotic solution at 
Mach number 1, but the part of the thickness and chord that Is more indic- 
ative of the portion of the airfoil that can influence the subsonic part 
of the flow field. This point has been discussed by Barish (ref. 27) who 
suggests the use of values of thickness t* and chord c* that are asso- 
ciated with the sonic point. This suggestion is based upon the observa- 
tion that numerical calculations indicate the asymptotic flow fields of 
wedges and certain cusped-nosed airfoils are nearly independent of the 
details of the shape, provided t* and c* are fixed. It follows, to 
whatever extent this observation is generally true, that the relation 
corresponding to equation (43 ) , but appropriate for nonaf finely, as well 
as affinely, related airfoils is of the following form: 



(44) 


where 


ACp* 


jX 2 (y+i )] 1/3 

t *2/ 3 


ACp 


H* = [^ 2 (7+1)t*] 1/3 H* 


h* = A 

c* 


T * = t* 

C* 


It is evident from inspection of sketches ( l) and (m) , as well as 
the more extensive set of results given in reference 2, that a substan- 
tial part of the discrepancies between the theoretical and exper imental 
results for Mach number 1 possesses properties consistent with those 
described above for the influence of the wind-tunnel walls. This corre- 
spondence prompts the presentation of sketch (p) which shows the varia- 
tion of ACp* with H* at free- stream Mach number 1 at the point on the 
airfoil surface at which sonic velocity occurs in free air. The v alu es 
for ACp* were computed using the theoretical values given in 
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Sketch (p) 


reference 2 for Cp and the 
experimental values_given in refer- 
ences 3 and 4 for C^ • The loca- 

*WT 

tion of the sonic point and hence 
the values of t* and c* are thus 
determined from the theoretical 
solution. Although the general 
procedure may "be open to some ques- 
tion and_scatter is inevitably large 
since ACp* is basically a small 

difference of larger quantities 
(some idea of the_uncertainty of 
the values for ACp* can be gained 


by consideration of the vertical line through each data point indicating 
the range of values associated with an arbitrary change of 1 percent in 
the local Mach number) , the results display a reasonably well-defined 
trend consistent with the ideas that ACp* does depend on H* and 


that | ACp* j diminishes with increasing H*. 


Pressure Drag 


Once the pressure distribution is known for a given airfoil, the 
pressure drag d can be calculated by use of the following relation; 


c d 



(45) 


It is convenient at this point to introduce a reduced drag coefficient c^ 
defined in terms of reduced quantities, thus 



(46) 


It is clear from the definitions of . Cp, t, Z, and 33 that the relation 
between c^ and c^ is 





1/3 


C d 


T 5/3 
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The variation of with g M has been computed for the symmetrical 
circular-arc airfoil. The results are summarized below and illustrated 
graphically in sketch (q) . Three classes of results are included in this 



Sketch (q) 


sketch and are distinguished by the type of line. The solid portion of 
the curve represents the values for Mach numbers up to that at which the 
shock wave has moved to 0.9 chord and is determi n ed by direct application 
of the pressure-distribution results for g < -0.985 given in table I 
and summarized in figure 3- The continuation of the curve to a value 
for g^ of -0.838 is determ in ed in a similar manner, but is indicated 
by a dashed line because the calculated pressure distribution is inaccu- 
rate in the vicinity of the trailing edge when the shock wave is at the 
trailing edge. The portion of the curve marked by long dashes represents 
the value 4.77 given in reference 2 for Mach numbers near unity. In all 
cases, the integrations required to determine the drag were evaluated 
using Simpson’s rule. Very fine intervals were used near the leading 
and trailing edges and the contributions of the regions in the immediate 
vicinity of the leading- and trailing-edge singularities were evaluated 
analytically in order to achieve the desired degree of accuracy. 

It can be seen by examination of the results shown in sketch (q) 
that c d is zero for g^ < -1.42 (corresponding to the critical Mach num- 
ber). It remains nearly zero as g M increases to -1.35* and then rises 
rapidly with further increases of g . The increase in drag is associated 
principally with the rearward movement of the shock wave across the chord 
and terminates abruptly when the shock reaches the trailing edge at 
g^ = -0.838. Although no values_for drag are given, it is evident from 
inspection of the variation of Cp with g ra indicated in figure 8 that 
further increases in g^ result In only small changes in e d , and that 
finally, at a value of g OT somewhat less than zero, ~ d becomes 
invariant with further changes in g . 
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Airfoils defined by 

Eq. CO 

Eq. (3) 


LJ 

L_"i 

" 1 

> 1 . 

f 1 


Point of maximum thickness, (xyfcL _ , . , n 

' Zmax ficient reliability to provide a clear 

. evaluation of the accuracy of the eal- 

Sketch (r) culated results shown in sketch (q). 

The most extensive results are those given by Michel, Marchaud, and 
Le Gallo in reference 3 for four airfoils of different thickness ratio. 
These results were obtained by integration of the measured pressure dis- 
tributions and suffer in quality, not only as a result of the factors 
discussed in the preceding section of this report, but also because of 
lack of sufficient pressure-distribution data for points near the leading 
and trailing edges. The results of Michel, Marchaud, and Le Gallo are 
presented, nevertheless, in sketch ( s ) together with the theoretical 
results from sketch (q). Also included in sketch (s) are the theoretical 
values for c^ for purely supersonic flow. The latter results are from 
reference 2 and were calculated by integration of pressure distributions 
obtained by application of a simplified form of simple wave theory that 
is consistent with the equations of transonic flow theory. 


The values for the drag of the 
other airfoils considered in the pres- 
_ ent investigation have been calculated 
only for Mach numbers near unity. These 
_ results, given originally in reference 2, 
are summarized in sketch (r) together 
_ with the corresponding results for 

related airfoils defined by equations (l) 
_ and (3) upon insertion of values for n 

_ . ranging from 1.50 to 6.05. 

6 

There are no experimental results 
for the pressure drag of symmetrical 
circular-arc airfoils that are of suf- 
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Sketch (s) 
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Examination of the results shown in sketch (s) reveals that the 
theoretical and experimental variations of (T d with exhibit the same 
general trends but are quantitatively quite different. Part of the dif- 
ference is undoubtedly associated with boundary- layer shock-wave inter- 
action. As discussed at greater length in reference 1, this phenomenon 
can have considerable influence on the drag, but is disregarded completely 
in the basic formulation of transonic flow theory. From the standpoint 
of drag, one of the principal effects of shock-wave bound ary- layer inter- 
action appears to be that the shock wave is moved to a more forward loca- 
tion along the chord. The pressures on the rear part of the airfoil are 
increased as a result and a substantial reduction in drag follows. These 
effects are largest at Mach numbers somewhat greater than the critical, 
but persist to a lesser degree for sonic and even purely supersonic flows. 
Because this phenomenon depends on Reynolds number and may be of greatly 
diminished importance at full-scale conditions, Michel, Marchaud, and 
Le Gallo introduced in the discussion of their experimental results for 
free-stream Mach numbers equal to or greater than unity the concept of 
"extrapolated drag" to represent the drag that would occur in the absenoe 
of separation. This quantity is calculated by consideration of a pressure 
distribution that differs from the experimental pressure distribution in 
the vicinity of the trailing edge as a result of the replacement of the 
pressures actually measured by those obtained by extrapolation of the 
trends indicated at stations upstream of the separation point. Accord- 
ingly, the values for extrapolated drag given by Michel, Marchaud, and 
Le G all o are also shown in sketch (s). As might be expected, the theo- 
retical values for drag are in better agreement with the values for extrap- 
olated drag than with those obtained directly from the actual measurements . 


Another part of the difference 
between the theoretical and experi- 
mental values for drag is associated 
with the higher critical Mach num- 
bers displayed by the experimental 
results in sketch (s) and is prob- 
ably to be attributed to the influ- 
ence of the wind-tunnel walls. Some 
idea of the magnitude of these 
effects can be gained by examination 
of sketch (t) , which shows the vari- 
ation of c d with for a symmet- 
rical circular-arc airfoil in an 
unbounded flow (or the variation 
of a, with for infinite H) 

<% °°WT 

together with the corresponding values for 



Sketch (t) 


WT 


and | 


°°WT 


computed by 


use of equations ( 39 ) and (40) for the same airfoil mounted in either a 
closed wind tunnel or open jet of such dimensions that H equals either 
or */2. These values of H are the same as considered in sketch (o) . 
Sketch (t) illustrates that the effect of the finite dime nsions of the 


test section is such that the values for c^ 


*WT 


remain equal to the 
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corresponding values for an airfoil in an unbounded flow, although 
experienced at different values of | , or Mach number. As a result, 
significant effects of the walls are observed in the variation of 

with e at Mach numbers somewhat greater than the critical. All such 
“WT 

differences disappear at Mach numbers near unity in an open jet or near 
choking in a solid-wall wind tunnel, however, since c^ is independent 
of | in this range. Examination of sketch (t) reveals again that the 
differences between the results for an unbounded flow and for an airfoil 
in an open jet are qualitatively similar to those that can be observed 
in sketch (s) between the theoretical results for an unbounded flow and 
the experimental results measured in a wind tunnel with slotted test 
section. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 22, 1957 
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TABLE I.- CALCULATED PRESSURE DISTRIBUTIONS 
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TABLE I . - CALCULATED PRESSURE DISTRIBUTIONS - Concluded 
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Figure 1.- Theoreti 
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Figure 2.~ Theoretical pressure distributions on an 

chord. 


airfoil with maximum thickness at 40-pereent 
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